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By using a density matrix renormalization group method, we investigate the ground-state properties 
of a one-dimensional three-orbital Hubbard model on the basis of a j-j coupling scheme. For B^^O, 
where B4 is a parameter to control cubic crystalline electric field effect, one orbital is itinerant, while 
other two are localized. Due to the competition between itinerant and localized natures, we obtain orbital 
ordering pattern which is sensitive to B^, leading to a characteristic change of Fsg quadrupole state into an 
incommensurate structure. At 52=0, all the three orbitals are degenerate, but we observe a peak at q=0 in 
Fsg quadrupole correlation, indicating a ferro-orbital state, and the peak at q=-K in dipole correlation, 
suggesting an antiferromagnetic state. We also discuss the effect of octupole on magnetic anisotropy. 
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In recent years, it has been widely recognized that the or- 
bital degree of freedom is one of the key ingredients to un- 
derstand the origin of diverse magnetic properties of strongly 
correlated d- and /-electron systems. ' Since in general, d- and 
/-orbitals are spatially anisotropic, orbital ordering leads to 
highly non-uniform spin-spin interactions under the orbital- 
ordered background, and there occur various kinds of mag- 
netic structures according to the type of orbital even in the 
same lattice structure. Here we note that for /-electron sys- 
tems, due care should be paid to treat spin and orbital degrees 
of freedom, since spin and orbital are tightly coupled with 
each other due to a strong spin-orbit interaction. In this situa- 
tion, total angular momentum provides a good quantum num- 
ber, instead of spin and orbital itself. Such a complex spin- 
orbital state can be classified from the viewpoint of multipole, 
which is represented by total angular momentum. 

Thus, in /-electron systems, ordering of multipole mo- 
ment has been discussed intensively both from experimen- 
tal and theoretical sides. Since among multipoles, dipole 
and quadrupole are corresponding to ordinary spin and or- 
bital in d-electron systems, respectively, ordering of dipole 
and/or quadrupole moments can be grasped by the d-electron- 
like theory. In this sense, octupole should be the first ex- 
otic multipole moment which we encounter in the solid state 
physics. In fact, a possibility of octupole ordering has been 
actively discussed for actual /-electron materials such as 
Ce^Lai.^Be,^"^ NpOa,^"'^ and SmRu4Pi2.''''^° 

In order to gain deep insight into multipole phenomena 
of /-electron systems beyond the phenomenological level, 
it is quite important to clarify multipole properties of a mi- 
croscopic model for /-electron materials. However, this is a 
rather difficult task due to the complexity of the /-electron 
model which should include charge, spin, and orbital degrees 
of freedom simultaneously. For the purpose, the construc- 
tion of a microscopic model for /-electron systems based on 
the j-j coupling scheme has been proposed.^' By analyzing 
such an /-electron model, for instance, the microscopic ori- 
gin of octupole ordering in Np02 has recently been estab- 
lished."'""^ Furthermore, magnetism and superconductivity of 
several kinds of /-electron materials have been studied from 
a microscopic viewpoint. 

However, multipole properties of the periodic systems have 



not been investigated satisfactorily at the microscopic level, 
even though their potential role has been emphasized in sev- 
eral kinds of /-electron materials. In the j-j coupling scheme, 
we can microscopically evaluate physical quantities related to 
multipoles by treating them as combined spin and orbital de- 
grees of freedom. In this paper, thus, we develop a numerical 
technique to study multipole properties from a microscopic 
viewpoint on the basis of the j-j coupling scheme. We report 
detailed microscopic analysis about multipole properties in 
a one-dimensional model for /-electron systems. It is found 
that incommensurability of quadrupole correlation function is 
changed by the effect of crystalline electric field (CEF). We 
also observe the change of magnetic anisotropy due to the ef- 
fect of octupole magnetic moment which belong to the same 
symmetry of dipole moment. Throughout this paper, we use 
such units as h=kB=l- 

First let us explain briefly the model construction for /- 
electron systems.^' In the j-j coupling scheme, we first in- 
clude the strong spin-orbit interaction, and consider only the 
lower sextet with j=5/2(=3 — 1/2), where j is the total angu- 
lar momentum. Taking the effects of CEF, Coulomb interac- 
tion, and /-electron hopping into consideration, we can ob- 
tain a multi-orbital Hubbard model for /-electron systems. 
In this paper, we study the ground-state properties of a one- 
dimensional system along the z direction with N sites includ- 
ing one electron per site under the cubic CEF effect. We note 
that due to the cibic symmetry, the results do not depend on 
the chain direction. 

The Hamiltonian is, then, given by 



(1/2) E i> 



(1) 



where fi^ is the annihilation operator for an / electron with 
the z component of the total angular momentum /i at site i. 
The hopping amplitudes are estimated from the overlap in- 
tegrals between /-electron wavefunctions, which are given 
by t^i^2 ±i/2~^^o with tQ={3/56){f fa), and zero for other 
cases. Here (//ct) is /-electron hopping amplitude through 
the sigma bond between nearest neighbor sites. Hereafter, we 
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(d) 



C'- 0.5 



(e) 



C'- 0.5 



take ^0=1 as the energy unit. 

In the cubic CEF term, -B^y is expressed by a cubic CEF 
parameter, B9. Note that due to the cubic CEF effect, the sex- 
tet with j=5/2 is spUt into doublet and Fg quartet states, ^ 
which are represented by 

|r7±) = yiTel ± 5/2) - 75761 T 3/2), (2) 

|rg±) = ± 5/2) + T 3/2), (3) 

|r|±) = I ± 1/2), (4) (C) 

where Fg and r| orbitals are introduced to distinguish two 
Kramers doublets in Fg, the subscript ± denotes pseudospin 
to label two states in each Kramers doublet, and in the 
right hand side is the eigenstate of the z component of the 
total angular momentum . The Coulomb integrals / in the j- 
j coupUng scheme are expressed by three Racah parameters, 
Eo, El, and E2.^^ 

We investigate the /-electron model (1) numerically by 
exploiting a density matrix renormaUzation group (DMRG) 
method,^' which has been developed for the numerical anal- 
ysis of low-dimensional strongly correlated d-electron sys- 
tems and quantum spin systems. In the present calculations, 
the finite-system algorithm is used for N=16 chains with the 
open boundary condition. The number of states kept for each 
block m is up to m=300, and the truncation error is estimated 
to be 10~^ at most. Note that the number of bases is 64 for 
one site due to the three orbitals, indicating that the size of the 
superblock Hilbert space becomes very large as x64^. To 
reduce the size of the Hilbert space, in general, it is useful to 
decompose the Hilbert space into a block-diagonal form by 
using symmetries of the Hamiltonian. We mention here that 
the Coulomb interaction is diagonal in terms of J^°', which 
is the total of jz, but the diagonal blocks classified by J^"' 
are mixed due to the hopping and CEF terms. In fact, we can 
make use of the integer 2J^°' (mod 8) as a good quantum 
number, as well as the total number of electrons. We empha- 
size that this is the first trial to apply the DMRG method to the 
microscopic /-electron model on the basis of the j-j coupling 
scheme. 

Let us first discuss the change of the orbital state due to the 
CEF effect. In Fig. 1(a), we show the dependence of the 
electron densities nr=J2i.cr{fIrafircr)/N in the Ff, r|, and 
T^ orbitals for Ea=5, Ei=2, and i?2=0.5, where fira- is the 
annihilation operator for an / electron with pseudospin a in 
orbital r at site i. At .64=0, all the three orbitals are degener- 
ate, but we find that the r| orbital is favorably occupied, while 
the r| and Fy orbitals are almost empty. This result can be in- 
tuitively understood from the free-electron energy dispersion, 
composed of a cosine band due to itinerant Fg orbitals and 
flat bands due to localized Fg and F7 orbitals, as shown in 
Fig. 1(b). At -64=0, where two flat bands are degenerate, the 
F| cosine band is occupied from the bottom up to the flat- 
band level, in the case of one electron per site. Namely, it is 
favorable to occupy itinerant Fg orbitals to gain kinetic en- 
ergy by electron hopping. For B^kO, electrons occupy the 
lower Fg level, and we find no remarkable change in the or- 
bital densities. On the other hand, for B2>0, electrons are 
forced to occupy the lower F7 level to gain CEF potential en- 
ergy. Note that in the free-electron energy dispersion, the Fy 
flat band goes down to the bottom of the F| cosine band at 
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Fig. 1. (a) Electron densities in tlie Fg, Fg, and F7 orbitals. The dotted 
lines denote the linear extrapolation to npy^l and n^b=0. (b) Schematic 
view of the free-electron energy dispersion. Local electron densities for (c) 
B°=0, (d) B$=0.01, and (e) S0=0.02. 



360S2=16io, i.e., B2/to=0.044. If we estimate \he value of 
B2/to that satisfies nrj=i from the present DMRG results, 
we obtain B2/to'^0.045, which is consistent with the estima- 
tion for the free-electron case. In order to obtain convergent 
results for B'l>Q.Q3, it is necessary to enlarge m and/or N in 
the DMRG calculations. Unfortunately, we have not yet fin- 
ished such DMRG calculations, mainly due to the limitation 
of CPU time and memory resources. In near future, we ob- 
tain DMRG results also for B2>0.03, which wiU be shown 
elsewhere. 

Within the present DMRG calculations, let us discuss the 
characteristics of the change in the orbital structure due to 
the CEF effect. For the purpose, we investigate the local elec- 
tron densities nir=J2^{fl^fira)- As shown in Fig. 1(c), r| 
orbitals are equally occupied in every site at ^4=0. On the 
other hand, when B^ is increased, we observe that electrons 
are accommodated in F7 orbitals in some sites, in which the 
electron density in the r7 orbital is practically one, as shown 
in Figs. 1(d) and 1(e). Moreover, such F7 sites are arranged 
at regular intervals, indroducing cuts of the chain due to the 
absence of electron hopping, while the F| sites form clusters 
to gain kinetic energy in the clusters. Note that the clusters of 
the Fg sites should weakly couple with each other through the 
Fy sites, since the F| orbital is also occupied in the F7 sites to 
some extent. Thus, the inconraiensurate orbital ordering pat- 
tern appears due to the competition between the itinerant Fg 
and localized F7 orbitals. 

Now we turn our attention to multipole properties to dis- 
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Tj multipole 



multipole operator 





dipole 


Jaux — Jx^ J4:Uy — '^y^ Jiuz — Jz 




quadrupole 


039„=(1/2)(2J2-J2-J2) 
039„=(^/2)(,/2-j2) 




quadrupole 


05gx=iV?.l2).JyJz 
05gy=(VS/2)JzJx 
05gz={V3/2)JxJy 


r2-[t 


octupole 


T2u={Vl5/6)JxJyJz 




octupole 


T4ux<m{2J^-Jx J^-Jpx) 

Tiuy=(M2){2jl-Jyj1-jljy) 

T4„, =( l/2)(2 J| - J, J2 - J2 Jz ) 




octupole 


nux=iVTE/6){JxJ^-JlJx) 






n,,y=(VTE/6){JyJl-JiJy) 

Tsuz=(VT5/6){JzJ^-Jpz) 







Table I. Multipole operators up to rank 3. The overline on the product de- 
notes the operation of taking all possible permutations in terms of cartesian 



components, e.g., JxJy=JxJy+JyJx- 
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Fig. 2. Osgy, quadrupole correlation functions for -64=0, 0.01, and 0.02. 
Inset shows the results in the region of small x in a magnified scale. 

CUSS the ground-state properties from the viewpoint of multi- 
pole. In order to clarify what types of multipole correlations 
develop, we evaluate the multipole correlation function 

Xr,(9) = (l/iV)^e^«('=-')(Xfcr,^ir,), (5) 

k,l 

where Xip is the multipole operator with the symbol X of 
multipole for the irreducible representation Tj in the cubic 
symmetry at site i, and (• • • ) denotes the average using the 
ground-state wavefunction. The definition of the multipole 
operators is given in Table I. Note that we consider the mul- 
tipole operators up to rank 3, and there are 15 types of mul- 
tipoles including three dipoles, five quadrupoles, and seven 
octupoles.^^ 

In Fig. 2, among 15 types of multipoles, we show a typical 
result of the O^gu quadrupole correlation function for several 
values of 5°. A sharp main peak at q=0 is observed, indi- 
cating a r| ferro-orbital state, and its amplitude is reduced 
with increasing B^, since the electron density in the r| or- 
bital is suppressed due to the CEF effect. On the other hand, 
as shown in the inset of Fig. 2, we observe a small peak at 
q=Tr for 3^=0, while the position of the peak moves away 
from g=7r to an incommensurate one with the increase of B^. 
This characteristic change in the incommensurability of the 
Fsg quadrupole state due to the CEF effect is in accordance 
with the change of the orbital structure. 

Let us here consider multipole properties also for the other 



Fig. 3. Multipole correlation functions of (a) dipoles, (b) F^g and Fsg 
quadrupoles, (c) T^u octupoles, and (d) T'2u and Fju octupoles at -84=0 
for Eo=5, El =2 and £2=0.5. The values at the peak q=iT of the correlation 
functions of the magnetic T^u (e) dipole and (f) octupole at 5^=0. The 
ratio Eq:Ei:E2 is fixed to 5:2:0.5 and the magnitude is varied. In the 
horizontal axis, Eq is normalized by the bandwidth VF=32to. 

components in more detail, with particular attention to the 
case of i?4=0. In Figs. 3(a)-(d), we show the multipole corre- 
lation functions. As shown in Fig. 3(a), we observe a peak at 
g=7r in each of the dipole correlation functions, suggest- 
ing an antiferromagnetic state. On the other hand, as shown 
in Fig. 3(b), the O^gu quadrupole correlation function has a 
peak at g=0 with large amplidute beyond the range of this 
figure (see also Fig. 2), indicating a T\ ferro-orbital state, as 
intensively discussed above. Note that the correlation func- 
tions of the quadrupole moments other than O^gu become flat 
without any significant structure, since they behave indepen- 
dently at each site. Note also that the correlation functions of 
O^gy and O^gz coincide with each other in spite of their dif- 
ferent symmetries, since in the atomic limit, the expectation 
values {Olgy) and {Olg^) taken by the r| state are found to 
be equivalent. In Figs. 3(c) and 3(d), we show octupole cor- 
relation functions. As observed in Fig. 3(c), the T^u octupole 
correlation functions have a peak at (7=77, in a way similar to 
the case for dipole, as observed in Fig. 3(a). In Fig. 3(d), 
we also find a peak at q='K in the T^ux and T^uy correlation 
functions. Note that the and T^uz correlation functions 
become flat and they agree well with each other, since again, 
(^2^ti)=(^5M^) in atomic Umit assuming that the ground 
state is Vi. 



4 J. Phys. Soc. Jpn. 

Here we remark the anisotropic behavior of the magnetic 
dipole and octupole moments, which belong to the same 
synnmetry. In Fig. 3(a), we find that the correlation function 
of Jiuz dipole moment is suppressed in comparison with 
those of J4ux and J4„y. This anisotropy is naturally under- 
stood from the fact that r|_|_ states are described by |±l/2) 
among j=5/2 multiplets.^^ Namely, dipole moment Ues in the 
x-y plane at a local level due to the Fg ferro-orbital state, and 
consequently, the correlation in the x-y plane exhibits a larger 
value than that along the z direction. On the other hand, as 
shown in Fig. 3(c), the correlation function of T4uz octupole 
moment is also reduced from those of T4ux and T4„y, which 
is the same behavior in the dipole. 

Let us also discuss the magnetic anisotropy in the strong- 
coupling hmit. In Figs. 3(e) and 3(f), we show the values 
at the peak q=TT of the r4u dipole and octupole correlation 
functions, respectively, as a function of the magnitude of the 
Coulomb interaction with the ratio Eq:Ei:E2 fixed to 5:2:0.5. 
It is observed in common that the correlation is enhanced due 
to the Coulomb interaction, but we find distinctive behavior 
in the anisotropy between dipole and octupole. In Fig. 3(e), 
we find that the J4uz dipole correlation remains smaller than 
those of J4ux and J4uy, even when the Coulomb interaction 
is increased, since the ground state is basically described by 
the r| ferro-orbital state still in the strong-coupling region. 
On the other hand, as shown in Fig. 3(f), the T4uz octupole 
correlation is much enhanced compared with those of T^ux 
and T4uy, and changes to be dominant in the strong-coupling 
region. Note that this anisotropy cannot be explained by anal- 
ogy of the anisotropy of the expectation values of octupole 
moments in the atomic limit, since {T4Tj_x)-{'^iuy)>{'^Iuz) 
for the Fg state. Namely, this anisotropy could be attributed 
to anisotropic interactions between F4„ octupole moments in 
the strong-coupling region.'^ Thus, the change of magnetic 
anisotropy appears in the strong-coupUng region due to the 
effect of the r4„ octupole moments. 

Finally, we note that with increasing B4, all the multipole 
correlations are suppressed. Since electrons in the F7 orbital 
caimot move in the chain, multipole moments fluctuate in- 
dependently at each site. On the other hand, we can observe 
incommensurate structures in the multipole correlation func- 
tions, as shown in Fig. 2. The detail of the dependence 
of the multipole correlations wiU be discussed elsewhere in 
future. 

In summary, we have studied multipole properties of /- 
electron systems on the basis of the three-orbital Hubbard 
model, by using the DMRG method. At .64=0, there appears 
the antiferromagnetic/ferro-orbital state, since itinerant T\ or- 
bitals are occupied to gain kinetic energy. With increasing B4, 
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electrons are acconmiodated in localized F7 orbitals, leading 
to the incommensurate orbital structure. We have also found 
the charactersitc behavior in magnetic anisotropy of the F4„ 
dipole and octupole. 
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and Technology of Japan. 
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